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Abstract 



We investigate the possibility of generating sizeable dipole radiations in relativistic 
theories of gravity. Optimal parameters to observe their effects through the orbital pe- 
riod decay of binary star systems are discussed. Constraints on gravitational couplings 



, beyond general relativity are derived 

1 Introduction 

G\ ■ 
O . 

General relativity predicts dominant quadrupole gravitational radiations. However, alter- 
natives to Einstein's theory actually extend the concept of gravitational charge beyond the 
Q" 1 mass, leading, in principle, to dipole radiations. 

& . ^ PaPCT ' ^ WiSh t0 inVeStigate P0SSibmty ° f g6nerating SiZeaWe diP ° le radi " 

tions, in the framework of two minimal, but conceptually important, extensions of general 

relativity. First, we emphasize that a violation of the equivalence principle provides us with 
a rather powerful indicator for gravitational dipole radiations. In section ||| we contemplate 
the dominance of dipole radiations in the Brans-Dicke scalar-tensor theory where the addi- 
tional gravitational charge is related to the gravitational binding energy of compact bodies. 
Section f| is devoted to a Kaluza-inspired vector-tensor theory where the new gravitational 
charge is proportional to the extra dimensional velocity. In the last section, we consider 
binary pulsars as indirect searching tools for gravitational dipole radiations. Monitoring 
the orbital period decrease rate of binary systems leads to limits on either the scalar or 
vector charges involved. In particular, we advocate realistic tight constraints through the 
analysis of existing neutron star - white dwarf binary pulsars. The resulting bounds on a 
scalar gravitational coupling might actually compete with the expected level of precision 
of future satellite experiments dedicated to probing this coupling. 
*FNRS Research Fellow (email: ywiaux@fyma.ucl.ac.be) 



1 



2 Gravitational charges and the equivalence principle 



In electromagnetism, the classical Larmor result ||, ||] for the dipole radiation rate due to 
cyclic motions is given by 

dt 6vrc 3 ' K > 

with D l = J2 a 1a r a,7 the dipole moment of the source and q, the electric charge. In Einstein's 
tensor theory, the analogous gravitational charge is the mass m itself. Consequently, the 
conservation of momentum in the non-relativistic limit implies that dipole radiations are 
forbidden in general relativity and the u 2 /c 2 -suppressed quadrupole radiations dominate, 
v being the average velocity in the (binary) system under scrutiny. The possibility to 
generate dipole radiations in relativistic theories of gravity beyond general relativity is 
therefore challenging from a phenomenological point of view. 

An intuitive way to ensure such a dipole radiation is to immerse Einstein's theory of 
gravity in a five dimensional spacetimeQ: a projection of the quadrupole motion on our 
4-dimensional spacetime corresponds then to a dipole one (see Fig.|l|). 




Figure 1: Dipole as projection of extra dimensional quadrupole. 
Indeed, the action for the geodesic motion of a test particle in five dimensions, 

8=-cj m ^, (2, 

reduces to 

f u 5 

S ~ - / {mc + (m— ) u^g 5fl } ds (3) 

in the weak field approximation, if ds = cdr is the 4-dimensional infinitesimal length 
(u M = dx^/dr and u 5 = dx 5 /dr) and for u 5 < c. The second term in this action corresponds 
to a vector gravitational interaction similar to electromagnetism. The appearance of a 
second gravitational charge (mu 5 /c) associated with the 4-dimensional vector field <75 M (x, t) 
circumvents the momentum conservation law. It allows for dipole radiations if particles 
x Our convention for the metric signature is as follows: (+, — , — , —...). 
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propagate in extra dimensions with constant intrinsic velocities. Such a possibility will be 
considered in details in section ||. For the time being, we simply notice that the reduced 
action implies the following equations of motion 

5 

mii a + mr;// ~ c 2 (ttA) {g^ a - g^ p )^ . (4) 
Consequently, the acceleration of a body in a freely falling frame is non zero since we obtain 

5 

a-#~ -c 2 (^)V#50, (5) 

in the non-relativistic limit. The right-hand side term being of gravitational nature, we 
face here a violation of the equivalence principle, i.e., a violation of the universality of free 
fall. We therefore have established a connection between the existence of dipole radiations 
and the violation of the equivalence principle, thus generalizing a conjecture made in ||. 

Let us take advantage of this link to single out another class of relevant metric theories, 
in which the mass of particles would explicitly depend on an additional scalar gravitational 
field 6: 



S = — c J m{4>)ds. (6) 
In the weak field approximation, this action for matter particles becomes 

S ~ — c j {m + (ms) <f)} ds , (7) 
with s = {dhxm/d(j))\Q, and the inferred equations of motion are 

mii a + ml^uV ^ c 2 (ms) ^{g m - u'V/c 2 ) . (8) 
In the non-relativistic limit, 

a — g ~ —c 2 sV4> , (9) 

and, once again, the existence of the second gravitational charge (ms) implies both a 
violation of the equivalence principle (see Eq.(||)) and dipole radiations (see Eq. ([?])). 

Additional scalar and vector gravitational fields seem to be unavoidable in quantum 
unifications of the four known fundamental interactions. In the next two sections, we will 
focus on popular prototypes of scalar-tensor and vector-tensor models of gravity resulting 
from String Theory and Supergravity, respectively. 



3 Brans-Dicke theories 

The Brans-Dicke theory [|[] is the minimal extension of general relativity including a second 
gravitational field. It is however already sophisticated enough, as we indicated in the 
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previous section, to encounter an equivalence principle violation and to allow for large 
dipole radiations. 



3.1 Scalar-tensor theory for compact bodies 

In the context of scalar-tensor theories, the mass of extended bodies will naturally depend 
on the scalar field's background value. Indeed, considering a freely falling frame (hence 
switching off the metric field), the scalar interaction survives. The internal structure of a 
body depends on it, and therefore its overall mass too. As proposed by Eardley JH, we wm 
treat a compact (i.e. self-gravitating) body in this framework as being poinlike, though 
with an explicit scalar dependence in its mass. The corresponding action (including now 
dynamics for the gravitational fields) of Brans-Dicke theory reads 

S bd = J ^gd 4 x - ^<S" a <S>, a ) - c I m($) ds , (10) 

where g is the absolute value of the determinant of the metric g^ u , $ is the massless 
scalar field and u, a dimensionless parameter. We are however only interested here in the 
linearized version of this action. The constant $o standing for the cosmological background 
value of the scalar field, let us define the metric and scalar perturbations as follows: 



$ , 

— - = l+aV2cK'ip. (11) 

The decoupling and canonical normalization^ of the corresponding kinetic terms in the 
resulting theory, 



^7 = ld ' x 



+ ( -d^ip 



(12) 



are implemented for k' = 87rG'/c 4 , (G' = l/$o being the bare tensor coupling constant) 
and 

a2 =2^T3' (13 » 
With such a normalization, the tensor and scalar propagators are given by 



fj.a u (3 i n/3 va _ p.v a/3 



q 2 + ie 



(14) 



2 Notice that the somewhat unusual additional y/c factor in the normalization of the gravitational fields 
allows us to transfer the whole dependence on the light velocity into the couplings. The 1/c counting will 
then be obvious when evaluating the intensity of radiation. 
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and 



A F (q) = ^—, (15) 
q l + ie 

respectively, for a transfer momentum q. Notice that the tensor propagator is the same as 
in general relativity. 

The linearized action for matter (assuming non-relativistic motion) gives the following 
interaction terms: 

where T^ v = pv^v u (p = mi 3 (x- x(t))) is the non-relativistic expression of the energy- 
momentum tensor for ordinary matter and T = pc 2 , its trace. The constant a is identified 
with the scalar coupling to matter and, consequently, we recover general relativity in the 
limit u) — » oo (see Eq.([0|)). The new gravitational charge s is the sensitivity of the variable 
mass with respect to the scalar field: 

_ dlnm 

S ^ <91n$|* ' { <) 

3.2 Strong equivalence principle violation 

The effective gravitational coupling G12 between two compact bodies (labeled by 1 and 2) 
here follows from the virtual exchange of a tensor (T) and a scalar (5*) gravitational fields. 
In the non-relativistic limit, the interaction potential is given by the scattering amplitude 
M = M T + M s of the two particles. In the limit of static bodies (T^ v = pc 2 5^5 u ), we 
get from Eqs.dH-H), 



M = Mc A [1 + a 2 (l - 2si)(l - 2s 2 )] G' 



m\m,2 



where N is a standard normalization factor, and q, the spatial component of the transfer 
momentum. The corresponding interaction potential is therefore the Newton potential 
with 

Gn = [l + a 2 (l - 2si)(l - 2s 2 )] G' . (19) 

In the case of test particles, this effective coupling constant reduces to 

G = (l + a 2 )G', (20) 

which defines the universal Newton coupling constant G in terms of the bare coupling G' . 
This relation allows us to interpret the charge s as the ratio of the internal gravitational 
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binding energy £l gr of the body to its mass m: 

<91nm 
S = ~dlnG 



= (21) 
mcr 

Before tackling the issue of dipole radiations, we would like to emphasize that, here we 
have in fact a violation of the strong equivalence principle, i.e., the violation of the universal 
geodesic motion for self-gravitating bodies only. Indeed, the free fall of a compact body 
towards a non compact source (si = 0) is slowed down proportionally to its own sensitivity 
(52): from Eqs,(|l~9|) and (|20|), we obtain 



G 12 = (l-£a 2 )G, (22) 

with 

the effective scalar coupling responsible for what is called the "Nordtvedt effect" |], ||]. This 
effect, anomalous with respect to general relativity, is responsible for a polarization towards 
the Sun of the Moon's orbit around the Earth. The Lunar Laser Ranging experiment, 
monitoring the motion of the Moon with an impressive current precision of the order of 
one centimeter, puts a tight bound on the amplitude of such a polarization. Working in 
the restrictive limit of the Brans-Dicke theory, this amounts to constraining the effective 
scalar coupling to extremely small values Q: 

£<610~ 4 . (24) 

3.3 Scalar dipole 



As already noticed in the first paragraph of this section, (see Eq.(16)), the source of the 
scalar field is m(l — 2s) rather than the mass m itself. Consequently, monopole radiations 
are still forbidden, but dipole ones (induced by the charge s) will appear in the same 
way as in electromagnetism. Let us therefore make the identification of the scalar Brans- 
Dicke interaction proportional to the sensitivity s with the vector interaction present in 
the Maxwell action 



S em = / <fx 



-F F^ v — A 



(25) 



4 ^ c 3/2' 

where j M = q5 3 (x — x(t)) v M . If we perform the substitution 

^- = Vl6vrG" as a (26) 
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in our basic Eq.(|]), taking into account a polarization correction factor 1/2 when identifying 
A^v^/c with (p, we readily obtain the scalar dipole radiation rate 

dEf p 



ir = ~3^ D U D U> (27) 

with D l , s = J d 3 x' psx' 1 . This is precisely the result of the more standard procedure of 
developping the radiative solution for cp in multipole moments. In the particular case of a 
binary source in Keplerian motion, the dipole moment is given by = miSifi +m,2S2^2 
and we get, to first significant order in the scalar coupling £ , the average rate 

1 f T ( dEg ip \ I dEf p \ 2G*n 2 M\ , 



T J Q \ dt \ dt 3 c 3 a 4 

,/u \ / \ Ik 

with \i and M, the reduced and total masses respectively, and a, the semi major axis of 
the classical orbit. The enhancement factor g(e) = (l + e 2 /2) / (l — e 2 ) 5 ^ 2 is a function of 
the orbital eccentricity e. 

The second kind of emission^ coming out in the context of these theories is the standard 
quadrupole tensor radiation |§, ||, || of general relativity (with again G' substituted for G): 



dE q T uad _ &_ 



-ryB Dij , (29) 



with D % 'j = J d^x' p {x n x n — ^5^r' 2 ). The Keplerian motion of a two-body system leads 
to^ the average rate 

(-^ 

with the enhancement factor /(e) = (l + 73/24e 2 + 37/96e 4 ) / (l - e 2 ) 7/2 . 

As already emphasized, dipole radiations (of order 0(l/c 3 )) could a priori dominate 
quadrupole ones (of order 0(l/c 5 )). But, before discussing phenomenological implications, 
let us turn our attention to the other previously mentioned theoretical frame where dipole 
radiations are also expected. 



4 Kaluza theories 

We have considered the possibility of dipole radiations arising from the coupling of a scalar 

gravitational field to the sensitivity of compact systems. Could this radiation rather come 

from the coupling of a vector field of gravitation to an extra dimensional velocity, as already 

suggested in section ^ ? 

3 We will neglect the scalar quadrupole (as well as a scalar monopole not even encountered by our 
linearized theory) as it would be of order 0(^/c 5 ), hence negligible compared to the tensor quadrupole, of 
order 0(l/c 5 ). 

4 Here G' and G12 are set equal to G. Once more, this corresponds to neglecting small corrections of 
order 0(£/c 5 ). 
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In 1921, Kaluza JO]] introduced for the first time the seminal concept of extra di- 
mensions in a superb attempt to unify gravitation and electromagnetism. In this theory, 
classical matter particles are allowed to propagate in a 5-dimensional spacetime and the 
momentum (mv 5 ) in the fifth direction is interpreted in terms of the electric charge (q) 



rn 



Vl6^G(- 



(31) 



However, this identification constrains the electric charge q to irrationally small values for 
elementary particles, which led to abandon temporarily the concept of extra dimensional 
spacetime. Here, we would like to preserve the gravitational nature of the induced vec- 
tor interaction and to analyze the observable consequences of the associated charge for 
macroscopic phenomena. 



4.1 Gravity in 4 + 5 dimensions 

The Einstein theory in 4 + 5 dimensions is given by the action 

1 



S, 



(4+5) 



2ck 



(4+<5) 



d^x^~gR + J d^x^~gC mat (4/, 



9M Nj 



(32) 



where k 



(4+5) = 87rG(4 +( 5)/c 4 . Here, Gu+s) stands for the (4 + <5)-dimensional coupling 
constant and §mn is the metric of the (4 + <5)-dimensional spacetime^]. If one defines the 
metric perturbation as c/mn = Vmn + ^/2cK(4+8)hMN , the corresponding linearized theory 
(which is the one of interest) reads 



nquad 
^(4+5) 



d M h AB d M h AB - -d M hd M h + -d M hd N h M N 



- L -d N h MN d R h MR 



— r n ±MN 

2c 



(33) 



where Tmn is the stress tensor for ordinary matter. For simplicity, we will assume that the 
extra dimensions are compactifled. The compactification on a torus (with identical radii R) 



implies that each field will decompose into modes h 



(n) , 



Only modes with the same n 



will couple to each other. Consequently, the integration over the 5 extra dimensions in the 
action simply leads to a volume factor V s which will renormalize the gravitational constant 
when going from 4 + 5 to 4 dimensions: G4 = G^+s/V 6 . The structure of the 4-dimensional 
interaction (fields and couplings) may be derived through the following "metric reduction" 
0: 



VMN + J^CK^ +S )h M N = VMN + V / 2~CK4 



A 



A 



fin 



-2aip mn 



(34) 



s Throughout this article, we will stick to the following indices conventions: capital indices run over all 
dimensions (0 < M,N... < 4 + 5), greek indices over the usual spacetime dimensions (0 < fi,v... < 3), 
latin jj-indices run over the usual spatial dimensions (1 < i,j... < 3), while latin mn-indices run over extra 
dimensions only (4 < m,n... < 4 + 8). 
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where K4 is the 4- dimensional coupling. The scalar cp = 5 mn ip mn is the trace of the 6(6+1) /2 
4-dimensional scalar fields, A^ n are 6 vector fields. The parameter a is again introduced to 
normalize the kinetic term for the scalar (p. This process of "metric reduction" is rather long 
and awkward, as one has to operate intricate field redefinitions before correctly identifying 
n = modes with massless mediators and modes with infinite (Kaluza-Klein) towers 

of massive mediators of the interaction. 

However, one may find out the field content (with correct couplings) much more easily 
via a "propagator reduction". Let us first consider how the (4 + <5)-interaction is mediated 
and decompose the propagator in a second step. In 4 + 6 dimensions, there is only one 
graviton. From Eq. (|33|) , it is straightforward to define its propagator as 

MA^NB 1 ^MB^NA _ _ 2 m MN m AB 
9(4H) + i£ 



G MN,AB (q) = V V +V^V V ^ (35) 



where q(4+s) is the momentum transfer in 4 + 6 dimensions. In order to understand the cou- 
plings to the 4-dimensional matter particles, we may once again consider the scattering of 
two particles by the gravitational interaction. For matter propagating in our 4-dimensional 
world, the corresponding amplitude reads 



,(2) 

L af3 

A/" £ G 4 T$G^(q)T®, (36) 

9(4+«) 



with 54, a given transfer momentum. From compactification, we know that each mode 
is associated with a mass: q 2 i+ ^ = q\ — m^c 2 , where m n c = '-jj-, The total scattering 
amplitude splits then into one massless and massive terms, 

M = Ml T + J2^n T , (37) 

n 

each of them corresponding to the interaction between matter fields through the propaga- 
tion of massless or massive gravitational mediators, respectively. 

Neglecting for a while the massive states, we may decompose the original propagator 
into massless tensor and scalar propagators in 4 dimensions, with fixed relative couplings: 



q\ + is 

+ 



G4 



T^^—Tv) . (38) 

<?4 T IS J 



6 + 2 

The corresponding effective action for gravitation in 4 dimensions is given by 

S^ T = S^ + S^, (39) 
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where the tensor part Sq is 



5, 



d 4 x 



2c M 



and the scalar one Sq reads 



(40) 



aS 
°0 



(fx 



2c V $ + 2 



(41) 



These are precisely the results obtained through the usual "metric reduction". They cor- 
respond to the linearized version of a Brans-Dicke scalar-tensor theory with strong scalar 
coupling, 



a 



S 1 

6+2 ~ 3 ' 



(42) 



if G4 = G' (see Eq.(|i6|) in the limit of test bodies (s = 0)). 

The massive modes[] are supposed to be heavy enough to ensure that the corresponding 
Yukawa contributions do not interfere with the experimentally verified 1/r law, above one 
millimiter. However, the resulting scalar-tensor theory is totally excluded by solar system 
observations. In fact, the well-known deflection of star light by the Sun is sufficient to 
rule out a strong scalar coupling. Indeed, the light bending by a heavy body is given 
by the scattering amplitude of a photon and a massive scalar, through the exchange of 
gravitational fields. The tracelessness of the photonic energy-momentum tensor (T™ = 
0) prevents from any coupling to the massless field (f (see Eq.(|38|)). Consequently, the 
predicted deflection angle @st is equal to the angle Qgr of general relativity, up to a 
sizeable rescaling of the effective coupling constant (see Eq.(pO|)): 



©ST = 77 &GR = .. 

G 1 + cr 



1 3 

&GR < T ®GR ■ 



(43) 



So, we have to invoke some mechanism by which the scalar ip also acquires a mass and 
decouples [jl^]. This amounts to sending the scalar coupling a to zero, providing therefore 
an obviously viable model. 

Such models with "gravity propagating in a (4 + 5)- dimensional bulk and ordinary 
6 For these modes, the "propagator reduction" method leads to: 



n \ n 



1 



q\ — m^c 2 + ie 



2 + S 



AT 



„Ca I//S 1 „fJ.l3„i^a 2 fit/ a/3 



-T 



q\ — m^c 2 + ie 



(2, 



with /9 2 = (2/3) [(5 — 1)/(S + 2)]. The corresponding 4-dimensional theory jll], |l2| encounters one massive 
graviton (whose tensor structure is known to be slightly different from the massless case given in 
and one massive scalar for each mode n. 
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matter confined on our 4-dimensional brane world" have recently been introduced to solve 



the hierarchy problem of elementary particle physics [lj, |ll], |12J. In the original models, 
the 6 extra dimensions are compactified, as we have chosen to do it here. Experimental 
constraints still allow these dimensions to be large, as long as they are smaller than one 



millimiter (see [14| and references therein). 

Back to our concern, looking for potential sources of dipole radiations, we abandon the 
so-called Large Extra Dimensional framework and allow for matter particles themselves to 
propagate also in the extra dimensions. Treating the extra dimensional velocity as a new 
gravitational charge (see Eq.(§)), we wonder whether its existence may lead to observable 
consequences in terms of dipole radiations. 

4.2 Equivalence principle violation 

For simplicity, we confine the motion of bodies to one extra direction, say x 5 . The unique 
vector interaction (V) easily comes out from the "propagator reduction" when considering 
once more the dominant scattering amplitude for two particles (1 and 2). Thanks to the 
extra dimensional coupling, this massless amplitude now reads 

M = M% + M% (44) 

where is the tensor part of M.q T given in Eq.(|38|), and 

= NcHGj^l ^^jf\ , (45) 
J ( v ^ qj + ie ^ a K ' 

where j?, = (pv 5 /c)v^ defines the vector current^ . In the non-relativistic limit, and 
considering static sources (T^ u = pc 2 <5q<5q), the scattering amplitude is then given by 

AWc 4 G(l-4#)^. (46) 
cr —q z 



It corresponds to a Newtonian interaction, with a renormalized gravitational constant 

C; 



Gi 2 = G(l-4^). (47) 



This is the mathematical expression of the equivalence principle violation in extra dimen- 
sional (vector-tensor) theories. Notice here some fundamental differences with the case of 
the Brans-Dicke (scalar-tensor) theory. On one hand, we have no clue of how the extra 
dimensional charge (v 5 /c) may a priori be linked to the nature of the body considered. In 
particular, there is no reason whatsoever to focus on compact bodies only. Eotvos-type 
experiments should also be considered in this context. They could tightly constrain the 

7 The conservation of this current readily follows from the linearized conservation equations in 4 + 5 
dimensions (T AIN N = 0), after integration over the whole compactified extra dimensional volume. 
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extra dimensional (vector) charge |15fl of laboratory-size objects. On the other hand, just 
as in the case of the strong equivalence principle violation in the Brans-Dicke theory, an 
anomalous acceleration is also expected for planetary-size bodies, though with two major 
conceptual differences. First, the two bodies (1 and 2) considered in this case need now 
to be charged (see Eq.(|47|)). And secondly, the additional interaction being of vector na- 
ture, this anomalous motion may either be slowed down or accelerated. On the contrary, 
the analytical expressions for scalar and vector dipole radiation rates turn out to be very 
similar. 

4.3 Vector dipole 

The effective vector-tensor action corresponding to the propagator reduction given in 
Eq.(|||) reads 

SX T = SZ + SX, (48) 



with 



S% = j d 4 x 



(49) 



Consequently, the dipole formula readily comes out from the identification of the extra 
dimensional velocity with the electromagnetic charge (see Eq. (|3l|) ) in our basic Eq.(||). 
The vector dipole radiation rate therefore reads 

dEy P _ 8 G ~ i ^ 

where the dipole moment is now D l = f d?x (pv 5 /c) x l . To first significant order in the 
velocities {G\i — ► G), the vector dipole radiations due to a Keplerian motion of a binary 
system is therefore given by 



\ dt 3 c 3 a 4 \ c c 

\ I K v 

where g(e) is the enhancement factor already defined for the scalar dipole (see Eq . (|28|) ) . 

5 Constraints from binary pulsars 

Binary pulsars are ideal laboratories for testing relativistic gravity, thanks to the extreme 
precision of their radio pulse. The most famous among them is the double neutron star 
system PS R B1913 + 16, discovered by Hulse and Taylor in 1974. The analysis of radio 
pulse arrival times provides us with a measurement of the decrease rate of the system's 
orbital period, in agreement with the general relativity prediction for gravitational radia- 



tions to within 0.4% [16]. Until now, the PSR £>1534 + 12 is the second and last example 



of a double neutron star system for which the measurement of the orbital period decrease 
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rate was found |17| to be in agreement with general relativity, though with a less impres- 
sive precisionF] of about 15%. However, in the last few years, there have been numerous 
discoveries |1, p(| ^l| of neutron star - white dwarf systems. Notice that the birth rate 
of these binaries is, according to the most recent models, at least as high as for double 
neutron star binaries. Prospects for relativistic parameter measurements in these systems, 
including orbital decay, are very encouraging, particularly in the case of PSR Jl 141 — 6545 
considered here below. 

Before coming up with more phenomenological considerations, we wish here to put 
forward the optimal conditions on the system's parameters under which gravitational dipole 
radiations dominate. The scalar (or vector) to tensor ratio of energy losses turns out 
to be simply expressed in terms of the model-dependent gravitational charges, and of 
the observable eccentricity and periastron shift (the latter including in fact the standard 
v 2 /c 2 -quadrupole suppression factor). From an astrophysical point of view, it is more 
convenient to introduce the corresponding ratio of the orbital period decrease rates induced 
by gravitational radiations. In the Keplerian approximation, these ratios are identical. For 



Brans-Dicke theories, Eqs. (|30|) and (|28|) , lead, to first order in £, to 



while for Kaluza-inspired theories, Eqs.(|30|) and (|5l|), lead, to first order in the vector 
charges, to 

[jSa) K \jSa) K 1 + T l" - 'I { A^J • (o3) 

The parameter 

AUJ = ^(1^) (54) 
stands for the general relativistic periastron precession in radians per revolution, while 

h(e) = 1 + e V 2 (55 ) 

1 j ~ 1 + 73/24e 2 + 37/96e 4 

is the reduction factor for eccentric orbits. The crucial parameter is therefore h(e)/ Au>: the 

more circular (h(e) — > 1) and the less relativistic the orbit (Au> — ► 0), the more important 

the dipole radiations (relative to the quadrupole ones). However, for maximum eccentricity 

(e = 1), the reduction factor h(e) is roughly 1/3 such that its impact on the question of 

dipole dominance is not dramatic. On the other hand, too small a relativistic precession 

8 More recently, a new test of relativistic gravity was achieved jl^] in the closest and brightest of the 
known binary pulsars in our galaxy, the neutron star - white dwarf PSR J0437 — 4715. However, this 
test is not provided by the analysis of radiation effects, but through the measurement of the (range and 
shape of) time delay undergone by radio pulses periodically passing close to the companion star. A similar 
measurement has also been achieved for PSRB1534:+ 12 leading to test the theory to within 1% jl^j , that 
is, a far better precision than what orbital decay analysis provides. 
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would correspond to a vanishing and therefore unobservable orbital decay. For further 
reference, the high eccentricity (e ~ 0.62) and relatively large periastron angular shift 
(Alo ~ 6.5 1CT 5 orw- 4.2 V -1 ) of PSR £1913 + 16 @ are not considered as optimal 
parameters for dipole dominance. 

Notice that these theoretical criteria selecting which binaries could give the best con- 
straints on dipole radiations are analogous to those defined for testing explicitly the strong 
equivalence principle in the strong field regime. It has been shown indeed |^2| that small- 
eccentricity long-orbital-period binary pulsars are also the best laboratories in which a 
polarized motion (towards the center of our galaxy) due to the strong equivalence princi- 
ple violation can be measured. However, in the light of the suggested alternative between 
scalar and vector gravitational interactions, any constraint resulting from such a direct 
analysis of the anomalous motion due to an equivalence principle violation is in principle 
ambiguous. Indeed, the weakness of the orbital polarization of a binary system (say, pulsar 
or Earth-Moon) might be attributed to an accidental cancellation between the scalar and 
vector perturbations rather than to the smallness of both the scalar and vector charges. 
Such an ambiguity is absent from an indirect analysis of dipole radiations since, in this case, 
the scalar and vector emission rates simply add. Hence, any bound on dipole radiations 
from binary pulsars provides a limit on (vf/c — v\jc) 2 and £(si — S2) 2 ■ 

5.1 Vector dipole in binary pulsars 

As already noticed, we don't know what really determines the vector charge (v 5 /c) of 
individual stars. However, its origin can be understood as follows. The extra dimensional 
momentum of elementary particles propagating in the compactified space is quantized. The 
allowed values for this "transverse" momentum actually define the masses of the associated 
Kaluza-Klein excitations in the effective 4-dimensional theory {q 2 A+ ^ = q\ — rr? n (?\ Hence, 
one might envisage the possibility that such excitations be trapped into stars through 
gravitational accretion, just as WIMP's ought to be. The vector charge for celestial bodies 
would then be defined by the fraction of the total mass m due to these Kaluza-Klein 
particles (t> 5 /c = ^ mxR/m)- But a dynamical model is needed in order to estimate such 
a fraction in neutron stars. 

5.2 Scalar dipole in neutron star - white dwarf binaries 

The sensitivity (s) of each body being related to its gravitational binding energy, the 
nature of the companion star is crucial for scalar dipole radiations. In particular, the 
fact that neutron stars have essentially identical sensitivities (s ns ~ 0.2) leaves little hope 
for discriminating between general relativity and Brans-Dicke with double neutron star 
binaries. Let us therefore turn to another class of binaries. 

Neutron star - white dwarf binaries are particularly interesting systems when searching 
for dipole radiation effects since a white dwarf is much less compact (s wc i ~ 10 -3 ) than 
a neutron star. Moreover, their orbital eccentricity is expected to be much smaller if the 
observed neutron star is the first-born star and was recycled by mass transfer from its 
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companion (TlJ. So, the determination of the orbital period decrease rate of a neutron star 
- white dwarf binary pulsar might strongly constrain the scalar coupling of gravitation. For 
illustration, we focus on the particular case of a recently discovered system, PSR J1141 — 
6545, which seems to be promising in terms of relativistic orbital parameter measurements. 
Given [19] its eccentricity (e ~ 0.17) and periastron angular shift (Aw ~ 5.0 10 -5 or 
lo ~ 5.3 °yr- v ), Eq.© implies 



(fS) -1 + 1.410^- (56) 

\ / |(1141-6545) 

An agreement to within 10% with the general relativity prediction^ would already give a 
constraint on the effective scalar coupling £, of the same order of magnitude as the best 
current limits given by solar system experiments. A more optimistic precision of 1% would 
provide us with the limit 

£ < 10~ 5 (w > 10 5 ) , (57) 

far better than the best current £ < 1.5 10~ 4 constraint and reaching the precision level 
expected from future satellite experiments dedicated to probing scalar gravitational cou- 
plings. Needless to say that an indirect measurement of dipole radiations would establish 
the existence of new gravitational charges. 



6 Conclusion 

Until now, the observed rates of relativistic orbital decay for binary pulsars are consistent 
with a pure quadrupole emission of gravitational waves, as predicted by general relativity. 
However, the appearance of dipole radiations is quite generic in effective models inferred 
from quantum theories that unify gravity with the other fundamental interactions. In a 
scalar-tensor theory of gravitation, a natural source for such radiations is the sensitivity (s) 
associated with compact bodies. We have shown how more precise measurements of the 
orbital parameters for peculiar binary pulsars consisting of a neutron star and a white dwarf 
would strongly constrain the scalar gravitational coupling (£) of the Brans-Dicke model. 
Nevertheless, if orbital decay rates ever suggest a non-vanishing anomalous coupling, we 
will never be able to disentangle the minimal scalar-tensor extension of Einstein's theory 
from a Kaluza-inspired vector-tensor one characterized by an extra dimensional charge 
(v 5 /c). The following interchange between the hypothetical gravitational charges 




9 The total mass of the system is M ~ 2.300 Mq. The mass function gives the bounds m p < 1.348 Mq 
and m c > 0.968 Mq for the mass of the pulsar (which actually corresponds to the statistical value for 
neutron star masses in known binary pulsars) and of the companion respectively. Those figures lead to a 
prediction for the orbital period decrease rate in general relativity of order T RG ~ —3.881 10~ 13 s/s. 
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leads indeed to identical dipole rates. Only a direct detection of gravitational waves could 
distinguish between longitudinal (scalar) and transverse (vector) polarizations, and would 
consequently fix the nature of the new charge involved. 
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